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ABSTRACT 

We present evidence for a metal-poor, [Fe/H]~ — L4 cr=0.2 dex, stellar halo component detectable 
at radii from lOkpc to 70kpc, in our nearest giant spiral neighbor, the Andromeda galaxy. This 
metal-poor sample underlies the recently-discovered extended rotating component, and has no de- 
tected metallicity gradient. This discovery uses a large sample of 9861 radial velocities of Red Giant 
Branch (RGB) stars obtained with the Keck-II telescope and DEIMOS spectrograph, with 827 stars 
with robust radial velocity measurements isolated kinematically to lie in the halo component primar- 
ily by windowing out the extended rotating component which dominates the photometric profile of 
Andromeda out to <50kpc (de-projected). The stars lie in 54 spectroscopic fields spread over an 8 
square degree region, and are expected to fairly sample the halo to a radius of '--^TOkpc. The halo 
sample shows no significant evidence for rotation. Fitting a simple model in which the velocity dis- 
persion of the component decreases with radius, we find a central velocity dispersion of 152kms~^ 
decreasing by — 0.90kms~^/kpc. By fitting a cosmologically-motivated NFW halo model to the halo 
stars we constrain the virial mass of M31 to be greater than 9.0 x 10^^ M© with 99% confidence. The 
properties of this halo component are very similar to that found in our Milky Way, revealing that 
these roughly equal mass galaxies may have led similar accretion and evolutionary paths in the early 
Universe. 

Subject headings: galaxies: spiral — galaxies: individual (Andromeda) — galaxies: individual (An- 
dromeda) — galaxies: evolution — Local Group 



1. INTRODUCTION 

One of the fundamental t enets of the Cold D ark 
Matter (CDM) paradigm ijWhite fc ReesI 11978(1 is 
that dark matter halos form hierarchically, via a 
series of mergers with smaller halos. This pre- 
diction gives rise to the natural expectation that 
the stellar hal o^ is formed from disrupted, accreted 
dwarf galaxies flnhnston Sis-urdsson fc HernnuistJII 9991: 
iHelmi & White' '1999a'; 'Bulloc k. Kravtsov fc Weinberd 
EOOl; Bullock fc Johnston 2004). 

CDM models predict that stellar halos form early, with 
the majority of the stars in the halo formed within a few 
relatively massive, 5 x 10^°Mq, dwarf irregular (dirr) - 
size dark matter halos. These were accreted '^10 Gyr 
in the past, likely prior to the formation of the bulge 
component via gaseous collapse, mergers and starbursts. 
Such models predict that the halo should be metal-poor 
and relatively smooth within the inner 50 kpc, where the 
tidal signatures of early stellar streams have been erased 
(^ullock fc Johnston 20055- The giant stellar streams in 
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the Milky Way ("Ibata et al."2001a';'Newberg et al.'l2003) 
and Andromeda (Ibata ct al. 2001b; McConnachic et 
|2004|) testify to the longevity of minor accretions which 
occurred since the metal poor halo formed. The average 
stellar halo density profile is expected to fall off with ra- 
dius more quickly than that of the dark matter because 
the stellar halo is formed from the most tightly bound 
material in inf ailing systems, while the majority of the 
accreted dark matter is stripped and deposited at larger 
radii. The discovery of fossil evidence of these accre- 
tion events via the identification of substructure within 
nearby stellar halos may provide the only direct evidence 
that structure form ation is hierarchical on small-scales 
(|Bullock fc Johnston 2005) . 

One of the most promising systems for the de- 
tailed decomposition and study of the stellar halo 
and its relation to the dark matter halo is the 
Andromed a galaxy, or M3 1 , ( Walterbos fc Kcnnicutd 
1988; Pritc het fc van den Bergh 19 8 8: Irwin ct a l. 200^ 
Geehan et al. 2006*: Icilbert et al.ll200^ . Often consid- 
ered the Milk y Way's twin, this galax y, at a distance 
of --785 kpc, iMcConnachie et al1l2005j) . hes sufficiently 
close that modern instrumentation can be used to mea- 
sure spectra of individual stars in the top few magnitudes 
of the red giant branch (RGB), providing access to their 
radial velocities and metallicities. In this paper we probe 
the outer stellar component of M31 using radial velocities 
and calcium II triplet (CaT) metallicities of RGB stars 
to study the overall kinematic characteristics. 

2. OBSERVATIONS AND SAMPLE DEFINITION 

Since 2002 we have undertaken a spectroscopic survey 
of RGB stars (selected by their optical colors and mag- 
nitudes) in M31 with the 10m W. M. Keck II tele scope 
and the DEIMOS spectrograph l|Faber et al.ll2003|) . The 
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54 spectroscopic fields shown in Figure 1 comprise 9861 
stars with measured radial velocities (shown in Figure 2), 
and represent our ongoing program to provide a com- 
plete kinematic and metallicity characterization of the 
disk and halo of M31, as well as the accreted fragments 
that now litter this environment. Details of the spectro- 
scopic observations and the selection of M31 RGB stars 
from the color-magnitude dia gram (CMP) hav e been 
presented in previous papers lllbata et alJ 120051 - 105; 
iMcConnachie et al. ''2003": Ferguson et al.ll2002 ^. Briefly, 
we used the standard DEIMOS multi-slit mode for low 
density outer fields, but adopted our own mini-slitlets 
approach targeting > 600 stars per mask in higher den- 
sity inner fields using small (1.5" long) slits, with data 
reduction using our own custom software pipeline (105). 
Stars are selected for masks based on the position of the 
red giant branch (RGB) in the color-magnitude diagram, 
with a limit of / < 22 for useful spectra. For most fields, 
except where a structure like And-IX ijGhanman et al.l 
120051) or the lHuxor et all l|200,'iD fuzzy clusters were tar- 
geted, stars were given equal priority from the RGB tip 
to 7=21.25, then lower weight for fainter objects. The 
S/N degrades to about 5 at the limit of useful data. A 
summary of all fields observed with DEIMOS, the radial 
velocity yields, and the breakdown of stars by kinemat- 
ics into halo, disk, and Galactic windows, is presented in 
Table 1. Representative individual spectra are shown in 
Figure 01 as a function of both RGB color and /-mag. 

Radial velocities measured from the calcium II triplet 
(CaT) can be used to kinematically isolate a likely sample 
of M31 halo stars. However, extracting a representative 
sample of M31 halo stars is hampered by the fact that 
many of the spectroscopic pointings have directly tar- 
geted the obvious pho tometric overdensities identified in 
iFerguson et alJ l)2002|) . Fortunately, our discovery (105) 
that an extended rotating component dominates the star 
counts in all these photometric substructures (with the 
exception of the g iant tidal stream to the South-East - 
l lbata et al.ll2001bl:fM cConn achie et al.ll2003Hlbata et all 
20041 iChao man et al. 2006: iGuhathakurta et al.ll2006al: 
iKalirai et al. 2006a), suggests a method to procure a 
sample of stars which is unbiased with respect to the 
photometric overdensities. A candidate halo component 
of stars can be extracted by windowing out in heliocentric 
velocity all stars which rotate with this extended compo- 
nent (shown in Figure 2). In other words, the photomet- 
ric representation of our windowed sample should be spa- 
tially smooth to first order (see 105 for further discussion 
of this point). We construct windows for the extended 
rotating component by adopting the average fit to the 
extended rotating component, (t„ = 50km/s Gaussian, 
from the kinematic model presented in 105, and consid- 
ering the window to span ±2cr (or ±100km/s in the disk- 
lag frame). ^ This is considered to be our "disk" sample, 
those stars exhibiting thin-disk- like kinematics. Between 
disk-lag velocities of —160 < viag < — lOOkms^^ we de- 
fine an additional component of intermediate kinematics, 
which we term "thick disk" in the subsequent discussion 
(to be discussed in detail in a subsequent paper). As 
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^ While the dispersion of the disk-hke component 
30km/s on average in individual fields (105), it is 



is about 
broader 

(~50km/s) in the coadded disk fields. This is partly a problem 
of the projection effects and the simplicity of the kinematic model 
currently adopted. 



Fig. 1. — The current coverage of our radial velocity survey 
of M31 with the DEIMOS spectrograph, in standard coordinates 
(x,f?). The outer ellipse shows a segment of a 55 kpc radius el- 
lipse flattened to c/a = 0.6, and the Major and minor axis are 
indicated with straight lines out to this ellipse. The inner ellipse 
corresponds to a disk of radius 2° (27 kpc), with the same incli- 
nation as the galactic disk. The distribution of 54 Keck/DEIMOS 
fields, comprising 9861 stars with measured radial velocities, are 
overlaid on the survey region, with field names tied to Table 1. 



the peak of the velocity distribution of the disk is ob- 
served to shift slightly between fields (typically by 20 
to 30kms~^) compared to the simple disk model of 105, 
the selection windows are adjusted accordingly by a small 
amount (hence the limits of the "disk" and "thick disk" 
selections in Figure 2 are slightly ragged). The disk- like 
population is revealed in Figure 13 as a Gaussian with 
disk-lag velocity of zero. Some skew in the distribution 
arises from the truncation in the North-East region of 
M31 imposed in order to reject contaminating stars from 
the Milky Way. This also affects the Vhei visualization 
of the disk in Figure giving rise to a much larger - 
470km/s disk peak (whereas the disk is actually sampled 
with roughly even numbers from our spectroscopy). 

The minor axis fields require an additional windowing 
criterion. For inner fields where the rotating components 
dominate, we have windowed out the "disk" component 
as it passes over v^e/=-300km/s (even though this re- 
moves any possible halo sample as well - clearly this is 
unavoidable). For outer fields where the halo dominates, 
we do not window at all. The first such field is F5 at 
20kpc on the SE minor axis (Figure 1), and all fields 
beyond that were considered to be "pure halo" . 

We do not explicitly distinguish a bulge component 
in this study, as it does not obviously distinguish itself 
in our predominantly major axis fields. We note that a 
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Fig. 2. — The radial velocity measurements from our 54 spectro- 
scopic fields are displayed as a function of distance along the Major 
axis, with the top panel showing heliocentric values, and the lower 
panel velocities corrected for the expected motion of the disk (using 
the model of 105). The foreground Milky Way stars (shown with 
small orange dots) are removed from our sample. The remaining 
stars are divided into a disk-like population (disk-corrected velocity 
close to zero; red dots), an intermediate thick-disk-like population 
(green dots), and a halo-like population (a broad component which 
is not rotating with the disk; black dots), and those metal-rich stars 
identified with the giant southern stream (blue dots). 



relatively metal-rich bulge-like cor nponent is reported a t 
minor axis radii from 12-20 kpc bv lKalirai et alJ l)2006bj) . 
consistent with other photometric studies of RGB stars 
at these radii (Durrell, Harris, & Pritchet 1994, 2001; 
Mould & Kristian 1986; Brown et al. 2003; Bellazzini et 
al. 2003; Worthey et al. 2005). However, this is tangen- 
tial to the present study. Very few of our fields lie within 
this radial range, and even fewer are used for measur- 
ing the halo spectroscopic metallicity (see §4 and Fig. 9 
which shows the fields where the numbers and quality of 
the spectra are sufficient for obtaining a halo metallicity 
measurement). As discussed above, our fir st minor axis 
field su itable for direct comparison with the lKalirai et al.l 
(j2006b) bulge is F5 at 20 kpc. This field has a halo 
sample metallicity of [Fe/H] = -1.18 (to b e presented in 
Fig. 9), and is consistent with the findings of lKalirai et al.l 
f2006b.) within errors of their bridging halo/bulge fields. 
To ensure that this result is not biased by Galactic dwarf 
star contamination, we have tested the field F5 with ve- 
locity cuts at both -300km/s and -170km/s, The result- 
ing measurements are [Fe/H] = -1.23 {v < — 300kms) 
and [Fe/H] = -1.18 {v < -170 kms). This is in the oppo- 
site sense from that expected if dwarf contaminants were 
significant (see below). Our current inability to clearly 
distinguish a possible bulge component in the inner-most 
regions of M31 may bias our measurements to slightly 
higher metallicity, although there is no statistically sig- 
nificant evidence that this is the case. 

The Galactic foreground (discussed below) and the Gi- 
ant Southern Stream are the last contaminants to win- 
dow from our stellar halo sample. The Giant Southern 
Stream has kinematics whi ch stand out as a h igh negative 
velocity peak in Figure 2 Ijibata et al.ll200^ . It appears 
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Fig. 3. — TOP: Histograms of the halo/disk populations shown 
in Figure |^ where thin and thick disk populations have been 
combined into a single histogram. Upper panel: The kine- 
matic sample selection, defined in the frame of disk-lag veloci- 
ties, with the disk-like population (shown in red) centered at zero. 
The halo-like population is chosen to be effectively non-rotating 
{^disk—lag > 160km/s). The halo sample is shown both with 
(dotted cyan) and without (solid blue) the metal-rich kinemati- 
cally coherent structure in the North-East region, which is likely 
to be associated with the giant southern stream (ChaDman et all 
120061) . Lower panel: The same samples are shown in the helio- 
centric frame, highlighting the broad halo-like Gaussian distribu- 
tion centered on the M31 systemic velocity. BOTTOM: Repre- 
sentative histograms of 9 individual fields spread throughout M31 
(labeled as in Fig. 0, and arranged roughly by location in M31. 
The double-peaked disk/Galaxy components in the top three pan- 
els (from the N-E region of M31 - F30, FIO, F8) reveal that our 
Vftei < —160 km/s Galactic windowing is indeed very conservative, 
but deemed necessary to procure a nearly pure halo sample. 



4 



in only seven of the fields, although another eight fields in 
the North-East region show a coherent kinematic struc- 
ture which is similarly metal-rich to the Giant Southern 
Stream, suggesting that this region represents a domi- 
nant component of the Stream's return on the other side 

of M31 (,Chapman et al 2006,^ . Nonetheless, the total 

number of stars in the Giant Stream is still relatively 
small, leaving the extended rotating component windows 
having the most significant effect on the halo sample. 

The sample selection in the disk-lag and heliocentric 
frames is shown in Figure |3 (where the black dots mark 
the halo stars of interest to this study) and Figure 01 
where the total distributions are shown. Figure |2I reveals 
the stellar halo windowed sample as a broad Gaussian 
distribution centered on Vhei '^-300km/s, along with the 
rotating thin and thick disk-like populations together as 
a single sample. Also shown in Figure Ol are representa- 
tive velocity histograms in nine individual fields spread 
throughout M31. The double-peaked disk/Galaxy com- 
ponents seen in the top three panels (from the N-E 
region of M31 - F30, FIO, F8) demonstrate that our 
v/ie/ < — 160km/s Galactic windowing is indeed very 
conservative, yet necessary to obtain an uncontaminated 
stellar halo sample. 

After all windowing procedures, 1207 stars remain in 
our halo sample. We then refine this sample to 827 stars 
by applying the following quality criteria (the same crite- 
ria as in 105): (i) The velocity of stars derived from two 
different sky-subtraction algorithms has to agree within 
20km/s : removes 35 stars (ii) Stars with continuum > 10 
counts (or about 40 e- total/pix) : removes 38 stars (Hi) 
Stars with ay < 20km/s : removes 138 stars (iv) Stars 
with cross-correlation peak < 0.05 : removes 189 stars. 
This sample can be considered to underlie the extended 
rotating structure of (105). 

2.1. Galactic dwarf contamination 

To complete our halo sampling, we must remove the 
foreground Galactic component (mostly dwarf stars), ob- 
served as a strong Gaussian peak centered on -61 km/s 
(105). The Galactic window is stationary (defined here as 
stars with Vhei > — 160kms^^), except for fields within 
the dense disk of M31 on the North-Eastern Major axis 
where we chose to remove only those stars with Vhei > 
Okms~^ as the source density is so high that Galac- 
tic contaminants are rare. Any Galactic contamination 
near the M31 systemic velocity (vhei = — 300kms~^) has 
to be Galactic halo, and the corresponding density of 
dwarf stars is minimal. The Be sanccon Galactic star- 
counts model ( Robin e t al.ll200'^ estimates 10 Galactic 
stars per M31 DEIMOS field of view (16.7' x 5') with 
Vhei < —160 km s^^) if all available candidates were ob- 
served. For Vhei < — 300kms~^ the Besangon model in- 
dicates the number of Galactic stars per DEIMOS field 
to be only 1.8 stars (typically <<1%). The Besangon 
prediction at velocities < -300 km/s is particularly rel- 
evant since more than 50% of our halo sample is cut 
at Yhei <-300km/s. In practice, we have typically ob- 
served <50% of the stars which are consistent with M31 
RGB colors available in the region covered by a DEIMOS 
mask. Given our choice of thousands of targets per 
field inside the inner ellipse in Figure 1 (and a maxi- 
mum multiplexing of 600 stars per DEIMOS mask 
with our mega-holes approach - see 105 for further dis- 



cussion), we estimate a negligible Galactic contamina- 
tion out to ~30kpc, and at most 10% contamination in 
our fields at larger radius. We also refer to the work of 
iGilbert et all 0006) which uses a multi-dimensional re- 
jection of Galactic dwarfs based on a variety of spectral 
features (most of which are unavailable to us, given the 
small extra cted wavelength regi me from our custom data 
reduction). IGilbert "eTall l|2006|) are therefore able to cal- 
ibrate the use of kinematic selection alone as a means 
to reject Galactic dwarfs, finding that <1% of stars cut 
at Vhei < — 300km/s will be dwarfs, in agreement with 
the Besangon model. Therefore our CMD selection of 
RGB stars coupled with kinematic sorting of the candi- 
date M31 stars by radial velocity is likely to remove the 
vast majority of Galactic contaminants. 

The windowing of stellar halo sample provides another 
check on Galactic contamination. Our North-Eastern 
stellar halo represents a more 'Galactic dwarf free' sam- 
ple as it is cut in velocity below -300 km/s, whereas 
the South- Western stellar halo is cut typically from - 
160km/s< Vhei <-300km/ s. As we might expect based 
on the Gi lbert et alJ l)2006j) work, we find almost no dif- 
ferences between our North-East and South- West hemi- 
sphere results, neither in kinematics (§3) nor in metallic- 
ity (§4), and thus it is clear that Galactic dwarfs are not 
significantly affecting our stellar halo results. 

A further reality check that our spectroscopically iden- 
tified stars are truly M31 RGB stars comes from the 
NaiA8i83,8i95 which is sensitive to surface gravity, and 
is accordingly undetectable in M31 RGB stars, but is as 
strong as the CaT absorption lines in Gal actic dwarfs 
with r elatively red colors (e.g., (Guhathakurta et "all 
l)2006aj) ^. We can check for this on a field by field ba- 
sis in combined spectra of stars, and where our spec- 
troscopic approach resulted in wavelength coverage to 

< 8200 A (the ~30% of fields where we employed the 
standard DEIMOS multi-slit masks). Our stacked spec- 
tra in all M31 components (extended disk, thick disk, 
halo, and Giant Stream) show very little evidence for 
NaiA8i83,8i95 absorption lines, whereas stars drawn ran- 
domly from -100 to zero velocity show strong Nai ab- 
sorption on average (Figure Isj . Direct measurement of 
the Nai equivalent widths in the halo and Galactic win- 
dows shows the contamination fraction is less than 15%, 
although this is of course a severe upper limit since there 
are large num bers of tru e M31 RGBs in the "Galactic" 
window ( Gilb ert" et alJ20 06'l diluting the Nai absorption, 
and weak Nai is also observed in RGBs.. This further 
corroborates the very low Galactic contamination in our 
M31 kincmatically-selected samples. 

3. RESULTS; KINEMATIC PROPERTIES OF THE STELLAR 

HALO 

To provide an overall kinematic characterization of the 
M31 halo stars, we first fitted a simple model with an 
isotropic velocity dispersion that is independent of ra- 
dius, but that is allowed to rotate as a solid body about 
the centre of the galaxy, which is assumed to have a 
systemic heliocentric velocity of — 300kms~^, and with 
the same axis of rotation as the disk. For this first fit, 
we manually removed 3 fields from the sample, due to 
the obvious presence of a kinematic halo substructure. 
These three contiguous fields lie on the North-Eastern 
major axis at the edge of the inner ellipse shown in Fig- 
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Fig. 4. — Representative spectra chosen from one of our survey fields (F3-w72), showing blue RGB stars with V — I = 1.3-1.8 (left) 
and red RGB stars with V — I = 1.8-2.5 (right), over the range of magnitudes for which we obtain good cross-correlations with the CaT 
template (/ < 22). In this field there are no RGB candidates brighter than / = 21, nor any good cross-correlations for stars fainter than 
/ = 22. 




Fig. 5. — Averaged spectra with wavelength coverage below 8200 A for typical DEIMOS fields. Left panels show halo samples from both 
North-East and South- West regions of MS, and therefore having velocity windows ranging from — 160km/s > v^^j > — 300km/s (F3,F4) 
to — 300km/s > Vftej > — 500km/s (F12,F9,F22). These fields show strong CaT absorption and minimal Nal>,8i83 3195 absorption. Similar 
results are found in all halo samples, providing a reality check on the minimal Galactic contamination to these samples at high negative 
velocities. The right panel shows the average of stars drawn from the v^el < — lOOkms"^ region thought to be highly contaminated by 
Galactic dwarf stars in throe of the same fields. These stars show a strong Nal absorption typical of Galactic dwarfs with redder colors 
(e.g., Guhathakurta et al. 2006), with weaker CaT absorption than seen in the metal-poor halo samples. There are between 10 and 16 
stars in each stack. 



ure 1, and the kinematic substructure is clearly visible on Figures 2 &|21with heliocentric velocity in the range 
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Fig. 6. — Likelihood contours (at Icr, 1u and 3cr intervals) 
of a simple rotating isothermal halo model, as a function of ve- 
locity dispersion and rotation. The asterix symbol denotes the 
position of the most likely parameters cjv = 126kms~^, £7 = 
— 0.08 km s^^/ kpc, where positive Q denotes motion in the same 
sense as the galactic disk. The data are clearly consistent with zero 
halo rotation at the Icr level. 



—400 < Vh < —350 km s ^ ( this d etection is is discussed 
further in iChanmaii et all f2006)). The model was fit 
simultaneously to the halo stars in all of the remaining 
51 spectroscopic field using a windowed maximum likeli- 
hood approach. The resulting distribution of likelihood 
as a function of the model dispersion and angular speed is 
shown in FigurelHI The best model parameters given this 
data are a halo velocity dispersion of cr„ = 126kms~^, 
and an angular velocity oifl = — 0.08kms^^/kpc. How- 
ever, the data are also consistent, at the la confidence 
level, with no rotation. 

Indeed, the halo shows no significant signs of rotation 
in any orientation. By dividing the spectroscopic fields 
into two sub-samples using a line intersecting the cen- 
ter of M31, we find the largest velocity differential as a 
function of azimuthal angle occurs near the Major axis 
where the NE half shows Vhei — —292 ± 8km/s from a 
best fit Gaussian to the average sample, while the SW 
half has v^e/ = — 303±6km/s. This insignificant rotation 
is in the same sen se as the extended rotating component 
ijibata et al.ll2005|) . and could for instance be a result of 
residual extended rotating component stars contaminat- 
ing our halo sample. For simplicity we therefore neglect 
rotation in the following, slightly more complex models. 

In Figure we show the result of fitting a model in 
which the velocity dispersion is allowed to vary linearly 
with projected radius R. The systemic velocity was again 
fixed at — 300km/s for stability due to a complex window- 
ing function. We proceeded by first calculating the veloc- 
ity dispersions in individual fields, given the adopted win- 
dow function, and assuming that the dispersion is con- 
stant and Gaussian at the position of the field. The re- 
sulting maximum-likelihood calculated values are shown 
in the upper panel of Figure H A single model with a 
linear dependence of velocity dispersion with radius was 
then fit to the full dataset. The measurements for in- 




^,(0) (km/s) 

Fig. 7. — Upper panel: Explicit single dispersion fits to the 
halo stars in the DEIMOS fields as a function of projected radius 
_R. The vertical bars show Itr uncertainties in these maximum- 
likelihood fits. The substantial variations in the fitted values and 
their uncertainties are to a large extent due to the low number of 
halo stars in most fields, and due to the different window functions 
in each field. Lower panel: Likelihood contours (at Itr, 2cr and 
3(T intervals) for a non-rotating halo model in which the velocity 
dispersion is allowed to decrease linearly with projected radius R. 
The most likely central velocity dispersion of this model is tT„ = 
152kms~^ which decreases with radius R as — 0.90 kms~^/ kpc 
(the straight line corresponding to these values is displayed in the 
upper panel). The fields in the upper panel that have a likelihood 
of less than 1% of being consistent with this line (displayed in red) 
were rejected in an iterative manner to produce this fit. 



dividual fields are not used in the analysis, presented 
only to clarify what the maximum-likelihood algorithm 
is latching onto. Proceeding in an iterative manner, we 
rejected those fields where the likelihood of the velocity 
dispersion conforming to the linear fit was less than 1%, 
and re-computed the fit. The likelihood contours of this 
model are shown in the lower panel of Figure[71 the most 
likely model has 

R 

aJR) = 152 - 0.90-— kms^^kpc^^ 
1 kpc 

At a distance of 10 kpc, the inferred velocity dispersion 
of the halo is therefore ~ 143 km s^^, while at a distance 
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of 60 kpc the inferred velocity dispersion has changed to 
~ 98kms~^. These results stand in contrast to the es- 
sentially is othermal halo measured for the inner Galaxy 
ijChiba fc Beers 2000); for M31, an isothermal model can 
be rejected at the ^ 3cr level. 

The individual field-by-field fits reveal that the over- 
all fit when we allow for a gradient in the dispersion is 
obviously dominated by the stars inside 40 kpc. While 
stars likely belonging to the Giant Stream have been re- 
moved from the halo sample, several stream fields stand 
out very clearly as not belonging to the trend by hav- 
ing a very low velocity d ispersion, as do fields in the NE 
shelf region (described in lFerguson et alJl2002l which also 
shows stars with giant stream kinematic structure) . The 
fact that these fields still stand out from an otherwise 
well-behaved trend suggests that our windowing of the 
giant stream is not extensive enough. However, at the 
present time we lack a full understanding of the distri- 
bution of the stream debris over the face of M31, forcing 
us to adopt this pragmatic approach. 

The above two models are essentially just phenomeno- 
logical descriptions of the data, so in the following we 
attempt a more physically-motivated model. Modern 
cosmology finds that th e dark matter halos of galaxies 
have a universal shape ijNavarro. Frenk fc Whit3IT997l 
hereafter NFW). Simulations show that halo stars are 
the remnants of ancient accreted proto-galaxies, but ones 
that were relatively massive. As such the stars deposited 
by these accretions are biassed towards being closer to 
the centre of the halo, and consequently the halo is less 
extensive than the dark matter halo ( Abadi ct al. 2006|). 
Nevertheless, the halo stars provide a reasonable first 
approximation to the dark matter, and we will therefore 
attempt to fit the population with an NFW model. This 
comparison is further motivated by the finding that the 
velocity dispersion of the halo sample decreases with ra- 
dius, which is expected from the NFW model. 

The NFW model is fully described by concentra- 
tion and virial mass parameters. Though the con- 
centration is on average dependent on the virial mass 
(pBuUock. Kravtsov & Weinberg 200f), there is a large 
scatter in this relation, and it is therefore safer to 
consider the full parameter space. We integrated the 
N FW model numer i cally, using the relations described 
in lLokas fc MamonI l)2001|) . to obtain the projected line 
of sight velocity dispersion for each field in our sur- 
vey, under the assumption that the velocity dispersion is 
isotropic. Note that the relative densities between fields 
were purposefully not specified, due to the fact our survey 
does not provide a complete census of stars in each field. 
We further adopted Hq — 70 km s~^ and a virial overden- 
sity parameter of 200. Figure |S1 shows the resulting likeli- 
hood contours (after rejection of the low-dispersion fields 
marked red in the upper panel of FigureCJ . Although the 
most likely model parameters are My = 9.8 x 10^"'^ Mq 
and c = 38, the likelihood contours are very extended in 
the concentration direction. The main reason for this is 
that our maximum-likelihood analysis method, as men- 
tioned before, does not have access to the information on 
the field-to-field density variations, whereas the concen- 
tration by its very nature relates to the radial variation 
in density. What density sensitivity we do have is due to 
the integration along the line of sight in each field. Cos- 
mological simulations indicate that concentration values 




6xl0" 8xl0" lo'^ 1.2xl0'^ 1.4xl0" 1.6xlo'^ l.Bxlo'^ 2xlo'^ 
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Fig. 8. — Likelihood contours (at Icr, 2cr and 3cr intervals) of an 
NFW model fit to our halo star sample, after rejection of the low- 
dispersion fields identified from the analysis presented in Figure 6. 
Here the most likely model has a virial mass of Wly = 9.8 X 10^^ Mq 
and concentration c = 38. Notice however, that the constraint on 
the model concentration provided by our analysis is not strong, 
and substantially lower concentration values are allowed (lower 
than c = 18 within the Icr boundary). The dashed line shows 
the c = 21 upper limit to the concentration of an M31-like galaxy 
argued by KlvDin. Zhao & Sommervillc (2002). Using this addi- 
tional constraint, the virial mass of the M31 halo must be greater 
than Mv = 9.0 X lO" Mq with 99% confidence. 



of c > 21 are v ery unlikely for t he dark halos of galax- 
ies such as M31 (|KlvDin. Zhao fc Sommervillc 2002); this 
limit is shown by a dashed line in Figurc|S| Imposing this 
additional constraint implies that the virial mass of the 
M31 halo must exceed My = 9.0 x 10" Mq with 99% 
confidence. This value for the halo mass agrees remark- 
ably well with the value calcula ted from the kinem atics of 
the giant stream around M31 (jibata et al.ll200j) . where 
we found a mass interior to 125 kpc (the apocenter of 
the stream) of M{R < 125 kpc) = 7.5^11 x 10" Mq. 
In contrast, the present constraint translates to M{R < 
125 kpc) > 7.2 X 10" Mq. 

4. RESULTS: METALLICITY PROPERTIES OF THE 
STELLAR HALO 

The photometrically dominant extended rotating com- 
ponent appeared to have a relatively high metallic- 
ity, [Fe/H]=-0.9 cr=0.2 dex (105), consistent with our 
larger sample of disk stars reported here, [Fe/H]=- 
1.0 (7=0.3 dex (where the error bars represent 
field to field dispersion in the measurements), ex- 
plaining why so man y disk-like stars have been 
found at large radius fReitzel fc G uhathakurtal I2002t 
iReitzel. Gu hathakurta fc Rich 2004). While attempts 
to model a more metal-rich halo have had some de- 
gree of succ ess under specialized evolutionary conditions 
l|Font et alJl200fiD . our CDM understanding of how the 
halo of a typical L* spiral galaxy forms suggest there 
must be a metal-poor relic of the early halo. 

But is there perhaps a metal-poor halo lurking un- 
derneath the relatively metsd-rich rotating component? 
Our first goal is to assess the metallicity of our velocity- 
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windowed sample of stars, to see if it shows properties 
distinct from the prominent extended rotating compo- 
nent which vastly outnumbers t he halo (on averag e by a 
factor ~ 6 X ) at radii < 50 kpc ijirwin et al. I l2005li. 

The DEIMOS spectra allow a measurement of the 
metallicity of the targeted stars from the equivalent 
widths (EW) of the Ca II triplet absorption lines, free 
from the model-dependency of photometric analysis. 
While the noise in individual spectra is typically too high 
to yield a useful measurement of the Ca II equivalent 
widths (whereby individual spectra can have vastly dif- 
ferent systematic errors due to the inopportune location 
of sky lines for certain heliocentric velocities), we esti- 
mate the average metallicity by stacking the RGB star 
spectra in individual fields. We shift each spectrum to 
zero velocity before stacking, normalizing the continuum 
to one using the inter-CaT continuum regions to fit a 
spline function, and finally weighting individual spectra 
by the continuum S/N (Figures & HHI) ■ Metallicities 
were measured in an identical way to 105, using the cor- 
relation of th e Call equivalent width with [Fe/H] as cal- 
ibrated by Ilutled£eiJiessei_&_&tetso (1997). 

For comparison, we derived photometric [Fe/H] values 
for the halo stars by interpolating between globular clus- 
ter fiducial sequences described in 105. The photometric 
properties of these metal-poor halo RGBs are entirely 
consistent with the spectroscopic properties, as shown in 
Figureim (we defer the model/age-dependent analyses of 
the stellar photometry to a foUowup paper). The largest 
variations occur in the N-E fields where the metal-rich 
giant stream structure is present at the halo window ve- 
locities (vhei < — 300km/s), which appear to be more 
metal-rich in the spectroscopic measurements. This may 
be indicative of the different ages/properties of the giant 
stream stars from the stellar halo. 

The average CaT triplet metallicity in stacked spectra 
with at least 10 high quality halo stars, is [Fe/H]=-1.4 
a = 0.2 (where the error quoted is simply the dispersion 
in the individual field measurements), much poorer than 
we find for the extended rotating component defined in 
the same manner [Fe/H]=-1.0 a = 0.3. There are five 
halo fields in the North-East which appear significantly 
more metal-rich than average and are not typical for the 
halo sample. These five fields all possess a kinemati- 
cally coherent sub-structure at Yhei ~ — 380km/s and 
with [Fe/H] '^-0. 7 which contaminates the halo-selected 
sample to various degrees. Including these sub-structure 
stars in the halo sample only marginally changes the re- 
suh ([Fe/H]=-1.36 a = 0.32 with, versus [Fe/H]=-1.40 
a = 0.21 without). This sub-structure is likely to be 
assoc iated with the giant southern stream (Ibat a~, al.l 
Upofb, 200^. 

Aside from these giant stream fields in the North-East, 
the metallicity measurements in each halo field are very 
similar, with no apparent radial gradient (FigurelHl). The 
three most outlying metallicities found are for fields D14 
and D12 (both with halo [Fe/H]— 1.7) and for F24 (with 
halo [Fe/H]=-0.9), none of which can obviously be at- 
tributed to selection effects or contamination. Fields D14 
and D12 are in the N-E hemisphere of M31, where as 
we've noted the halo sample cut with Vkei < — 300km/s 
results in < 1% Galactic contamination. While a small 
fraction of Galactic contamination is expected for F24 in 
the S-E hemisphere, this would in fact make the spec- 
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Fig. 9. — The metallicity measurements from stacked spectra 
are shown for the halo sample (left panels) and disk sample (right 
panel). The upper panels show the spatial location of the fields, 
with the ellipses denoting the same limits as in Figure 1. The 
smaller number of fields compared to Figure 1 is due to the require- 
ment that at least ten well-measured spectra of the same galactic 
component be present in a field to provide a metallicity measure- 
ment. The five metal-rich halo fields in the North-East are not 
typical for the halo sample, and are clearly the result of contami- 
nation from a [Fe/H]~-0.7 dex kinematically coherent structure at 
380km/s which appears in all these fields, likely to be associated 



with the giant southern stream IChapman et 



Kel y to o e 
ani200a . 



troscopic [Fe/H] appear more poor, not rich. These 
fields may therefore host accreted substructures which 
are not as well phase mixed as the halo in general. Our 
findings therefore suggest that there is a metal-poor, 
non-rotating halo component in M31, with comparable 
metallicitv to th e Milky Way's non-rotating stellar halo 
(IChiba & Beers 2000. 2001). 

The accuracy of our M31 halo metallicity measurement 
is limited by the continuum fitting in these relatively low 
signal-to-noise spectra, with additional uncertainty from 
systematic effects such as residual sky. The large dis- 
tance of M31 implies a small uncertainty for the horizon- 
tal branch magnitude, even if the stars are drawn from 
relative extrema in the galaxy. Furthermore, the RMS 
dispersion in T^-band magnitudes of the total sample of 
halo stars is very small (<0.5 mag). These effects dimin- 
ish the impact on the derived metallicity of an uncer- 
tainty in distance in the RGB population we are observ- 
ing in our sample. A formal error for each fit is obtained 
from the signal-to-noise of the CaT lines, with an extra 
error estimated for the uncertainty in sky subtraction. 
The typical uncertainty is ^ 0.3 dex. 

A concern is that the average CaT metallicity in a 
given field may be dominated by a few outlying spec- 
tra, although this should only be a factor for increasing 
the metallicity (increasing the EW). This potential sys- 
tematic error can be tested by dividing up the fields into 
halves and comparing metallicity measurements. Gen- 
erally subsamples have metallicities consistent to within 
0.2 dex. However, a danger with this approach is that for 
a small enough sample, variations will be present due to 
the random sampling of the underlying metallicity dis- 
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Fig. 10. — Representative combined spectra (weighted by the S/N of individual spectra in the stack) for kinematically selected halo 
stars in 14 of the 23 fields shown in Fig.|5j along with the total stacked halo spectrum at the top exhibiting an average [Fe/H]=-1.4 a=0.2 
dex, the error bar indicating the field to field RMS variation (whereas measurement errors in the stacked spectra are much smaller). Field 
names and metallicities are indicated on the left and right respectively. The spectra have been normalized to a continuum of one with the 

o o 

Ordinate crossing at zero, and smoothed to the DEIMOS instrumental resolution. The three Ca II triplet lines at 8498.02 A, 8542.09 A, 

o 

and 8662.14 A are the most obvious features of these spectra, although many other real features are present. 



tribution. 

Our sample of RGB stars isolated to have halo kine- 
matics (and importantly windowing out the extended 
rotating component) have allowed a census of the halo 
along both Major and minor axis fields, and underly- 
ing the range of photometric substructures observed. 
The uniformly low ([Fe/H] = -1.4 a=0.2) metallic- 
ity we measure across the fields reflects a stellar halo 
component which has a distinctly different origin from 
the extended rotating component. Our discovery is 
complemented by the independent discovery of a halo 
component with metal poor ([Fc/H ] = —1.5) stars ly- 
ing at large radius (60-150 kpc) (,Kalirai et alJ 120060: 
IGilbert et al.l l2006). Their outer halo component and 
bridging halo/bulge component at smaller radii are sim- 
ilarly metal-poor to the halo we describe h ere. However, 
the 47 pure halo stars at i? > 60 kpc in iKalirai et al.l 



l)2006b(l and lGilbert et all l|2006|) lack suflicient statistics 
to ascertain whether they are consistent with an outer 
extension of the kinematically characterized, pressure- 
supported, non-rotating component we describe in this 
paper. 

The larger field to field halo metallicity d ispersion (a ~ 
0.3) f ound by IKalirai et all l|2006hD and IGilbert et all 
l)2006j) is suggestive of an increasingly inhomogeneous 
halo component with radius, and is consistent with the 
dramatic structural variations from wide-field photomet- 
ric surveys of M31's outer halo by Ibata et al. (2006) and 
iMartin et all l|2006il . These finding s are also consistent 
with th e much less phase-mixed out er halo in the simula- 
tions of iBuUock fc JohnstoiJ l)2005|) . suggesting that the 
outer halo may truly represent a component with a much 
longer assembly and relaxation time than the inner halo. 
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Fig. 11. — Comparison of the spectroscopic and photometric 
[Fc/H] derived for RGBs in the halo sample. Points are shown for 
all fields with at least 20 stars. The line denotes equality. Note that 
the photometric [Fe/H] is highly age dependent; the good agree- 
ment with the spectroscopic [Fe/H] suggests the adopted ~13Gyr 
ages adopted from the globular cluster fiducials is reasonable for 
the halo stars. The metal-rich points correspond to fields in the 
N-E region of the extended disk where the giant southern stream 
(Ibata et al. 2001) passes through at velocities within the halo 
window. 



5. DISCUSSION 

While our study has uncovered conclusive evidence 
for a pressure supported metal-poor halo in JVI31 within 
the inner 70 kpc (complemented by the co-discovery of 
a metal-poor halo in Kalirai et al. 2006 and Gilbert 
et al. 2006), evidence for met al-poor stars i n An - 
dromeda's halo reaches back to iHolland et al.l l)1996t) . 
who used deep imaging in the outer halo to probe the 
halo CIVID. Even earlier results from RR Lyrae studies 
suggested tha t old, metal-poor stars might b e present 
in the halo ((Pritchet fc van den Berghl I1988D . IVIetal- 
poor RGB stars were fi rst identified spectroscopically 
in the halo using Keck fReitzel fc Guha thakurtal l2002t 
^eitzcl. Guhathakurta & Rich 2004), with speculations 
that differences appeared to exist between subsets of the 
RGBs in their sample, finding a large range in metallic- 
ities, and suggesting that disk-like stars were found at 
radial distances as large as 34 kpc (de-projected). How- 
ever, with only two isolated pointings and small numbers 
of stars (23 and 29 respectively), these studies had no 
means to characterize the outer structural components 
in M31. 

While our Keck spectroscopic data do not suffer from 
age/metallicity degeneracy as do the earliest photomet- 
ric studies of the halo, the most important ingredient for 
our discovery is the detailed knowledge of the extended 
rotating component which allows for the separation of 



populations by their kinematics. Thus when stars are se- 
lected in a similar way to the halo star sel ection in Galac- 
tic st udies (i.e. by extreme kinematics - iChiba fc BeersI 
12000(1 . a similar [Fe/H]~-1.4 metallicity is found. 

Since JVI31 and the JVIilky Way reveal so many differ- 
ences in their evolutionary histories and environments 
(e.g.. Bulge size, thin disk scale, number and types 
of globular clusters, scale siz es and types of satel- 
lite dw arf galaxies - e.g., iWalterbos fc K ennicT ^ 
1988; Barmbv et al. ^20^ iHuxor et aP l200a 
McConnachie et al, 2006a Jj)7 it is of considerable 
interest that they seem to have developed similar 
metal-poor stellar components in their outskirts, likely 
near their for mation epoch. If models of halo for mation 
are correct (Bullock & Johnston 2005; Renda et alJ 
12005: Fo LLt ct al- 2,006) , this similarity in metallicity and 
dispersion, along wi t h the similar dark matt er masses 
l|Evans fc WilkinsonI l2000l llbata et all 1200^ . suggest 
that both galaxies must have originated in comparable 
overdensity contrasts and perhaps attracted similar 
numbers and masses of merger accretions early on. 

Our finding of a low halo metallicity in JVI31 con- 
trasts, however, with recent measurements of the stel- 
lar halo metallicity distri bution in a sample of m ore dis- 
tant edge-on galaxies bv lMouhcine et aLl l|2005() . These 
authors derive photometric metallicity estimates from 
HST/WFPC2 observations and find a strong correla- 
tion between halo [Fe/H] and galaxy luminosity (Fig- 
ure El- For a luminosity similar to 1VI31 or the JVIilky 
Way, their relation predicts a mean halo abundance of 
[Fe/H]'^ —0.8, significantly more metal-rich than that 
found here . It is importan t to po int out that the fields 
studied bv HVTouhcine et all l|2005D are assumed to probe 
the stellar halo on the basis of their projected radii alone 
but this assumption has not been verified by either struc- 
tural decomposition, or kinematics. We have shown that 
when halo stars are selected kinematically - arguably 
the best way in which to do so ~ then the metallicities 
inferred for the Milky Way and 1VI31 are both similarily 
low. When combined with the recent discovery of metal- 
poor, [Fc/H]^ —1.5, kinematically-s elected halo stars in 
the lo wer mass spiral galaxy 1VI33 ijlVIcConnachie et alJ 
l2006c|l . this suggests that the true underlying stellar halo 
metallicity in spiral galaxies may have very little depen- 
dence on galaxy luminosity. 

In Fig. El we place our kinematically-selected halo 
metall icity results for JVI31 and M33 in the distribution 
of Renda et al.l l(2005fl . Their results suggest that galax- 
ies with metal-rich stellar halos at z = have a longer 
formation history, whereas galaxies with a more metal- 
poor stellar halo at z = have a shorter assembly, con- 
sistent with their halo [0/Fe] predictions. By kinemati- 
cally selecting only the non-rotating halo components of 
the Local Group spirals, it is clearly a sho r ter as sembly 
time, primordial halo in the iRenda et al.l (12005*) simu- 
lations that would be isolated. However, ^^nda et a^ 
(2005) define the halo as all the stellar particles in the 
simulation at a projected distance R > 15 kpc. Clearly 
a direct application of this definition to 1VI31 would re- 
sult in an average [Fe/H]~-1, being drawn primarily from 
the dominant ext ended rotating comp onent (105). Note, 
however that the IRenda et al.l (|2005) simulations (Fig- 
ure EJ suggest that even defining the total halo in this 
way, finding a halo metallicity as rich as [Fe/H]~-1 is 
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rently insufficient to show any kinematic differences or 
even allow robust metallicity measurements for the red 
and blue RGBs. Future expansions of the dataset may 
aUow the characterization of this population. 

Finally, it is of interest to understand how our new 
metal-poor component connects to the extended 
halo componen t discovered p hotom etrically along the 
minor axis by llrwin et al. I l)2005t). and also discov- 
ered through spectroscopic confirmation of halo-like, 
metal-poor stars at very large radius (60 — 150kpc) 
ijKalirai et al. 2006bt [Gilbert et al.ll2006j) . Are these aU 
part of the same halo component, and is there a gra- 
dient in kinematics and metallicity? Again, it is diffi- 
cult to assess, although the fact that massive satellites 
can be more easily accreted on prograde orbits (con- 
taminating the disk, but leaving the spherical halo in- 
tact), suggests that such an extension of the halo would 
be physically motivated. Tentative evidence from this 
work and from IKalirai et aU l)2006bj) for a more hetero- 
geneous outer halo component distinct from the more 
phase mixed inner halo component is broadly consistent 
with the expectations from theory ijBuUock fc JohnstonI 



Fig. 12. — The stallar halo metallicity - total galactic V— band 
luminosity relation. Our kinematically selected metallicity mea- 
surements for M31 (this work) and for M33 QylcConnachie et aLl 
I2006crt are shown as revised points connected by lines to litera- 
ture values. Ope n triangles represent th e HST data on other spi- 
ral galaxies from IMouhcine et al] 120051) ; op en squares represen t 
additional literature points extracted from IHenda et al.l 120051) . 
Filled squares correspond to the peak of the metallicity distribu- 
tion function (MDF) in halo simulations of Rcnda et al. (2005), 
where 68% of stellar particle in the MDF are enclosed within 
±0.7dex. The ensembles at different luminosity ranges represent 
all the stellar particles in the simulation at a projected distance 
R > 15kpc, and corresponding to total masses (including dark 
matter) of lell,5ell, lel2,5el2M0 respectively. 



quite rare for virial masses as la rge as M31 IO-'^^Mq). 
By contrast, in the iFont et aU Ij2006^) simulations find- 
ing an early assembled stellar halo with a lO^^M© virial 
mass and an [Fe/H]~-1.4 (like the MW and M31) occurs 
in only '^10% of the accretion histories, suggesting that 
the MW and M31 halos are actually quite rare for L* 
spirals. 

However, our discovery should not be interpreted as 
claiming that all components of the R < 70kpc stel- 
lar halo are metal-poor, but simply that a non-rotating 
metal-poor component to M31 exists which is similar in 
[Fe/H] to the halo in the Galaxy. There may also be 
a rotating halo sub-component which is more metal-rich 
which could effectively be windowed out of our current 
sample. Contamination from extended rotating compo- 
nent stars both in velocity and metallicity makes a de- 
tailed compa rison difficult. iFerguson et al.l 1)2002(1 and 
llrwin et al. I (|^005) noted that there is very clearly a 
persistent population with red colors (suggesting stars 
which may be metal-rich or young) along the minor axis, 
and beyond the region that 105 have characterized as 
the extended rotating component. While measuring the 
spectroscopic properties of this population would be in- 
teresting, the statistics in our minor axis fields are cur- 



6. CONCLUSIONS 

We have presented evidence for a metal-poor, [Fe/H] = 
— 1.4, cr — 0.2 dex, stellar halo component in M31, by 
kinematically isolating 827 non-rotating stars from our 
survey of ~10,000 stars with radial velocities from Keck- 
II/DEIMOS. The halo component is detectable at radii 
from lOkpc to TOkpc, and underlies a more dominant 
extended rotating component. 

This metal-poor halo component has no detectable 
metallicity gradient, and is consistent with an early 
and rapid formation period as suggested by the simula- 
tions of Renda et al] l|2005j) . The more metal-rich stellar 
halo components observed with HST bv IMouhcine et all 
may be indicative of components accreted later 
over longer assembly times. 

Our analysis has also allowed us to constrain the kine- 
matic properties of the M31 halo. The windowed halo 
sample shows no evidence for rotation about any axis, 
though there is clearly a drop in the velocity dispersion 
of the population with distance away from the galaxy 
center. Fitting an NFW model to these data and using 
the cosmologically-motivated constraint that the model 
concentration c < 21, provides a lower limit to the virial 
mass of the halo of My > 9.0x10" M© (99% confidence) . 
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TABLE 1 

Log of Keck/DEIMOS spectroscopic observations in the Andromeda galaxy. 
Positions are in degrees away from 00h42m44.3s +41d16m09s, J2000.0. 



number field name Xki Eta no. stars" no. halo* no.disk'^ no. gal'' obs.run 



1 


Fl 


wll 


-1.398 


-1.670 


108 


17 


76 


15 


sept.2002 


2 


F2 


50Disk 


-0.968 


-1.682 


266 


26 


184 


56 


sept.2004 


3 


F3 


w72 


-1.463 


-0.989 


95 


6 


72 


17 


sept.2002 


4 


F4 


w80 


-1.120 


0.020 


79 


6 


49 


24 


sept.2002 


5 


F5 


57Halo 


1.087 


-0.919 


113 


6 


24 


83 


sept.2004 


6 


F6 


lODisk 


-0.443 


0.462 


226 


4 


146 


76 


sept.2004 


7 


F7 


59Halo 


2.364 


-1.488 


96 


3 


13 


80 


sept.2004 


8 


F8 


54Blob 


-0.532 


1.061 


206 


8 


113 


85 


sept.2004 


9 


F9 


w91 


-0.297 


1.288 


96 


10 


60 


26 


sept.2002 


10 


FIG 


56Blob 


1.398 


0.187 


191 


12 


97 


82 


sept.2004 


11 


Ell 


w42 


1.380 


0.395 


333 


85 


49 


199 


sept.2002 


12 


E12 


w95 


0.530 


1.713 


122 


20 


38 


64 


sept.2002 


13 


E13 


06Disk 


1.392 


1.453 


176 


57 


13 


106 


sept.2003 


14 


E14 


14Disk 


1.041 


1.807 


76 


9 


2 


65 


sept.2004 


15 


E15 


53Blob 


1.943 


1.997 


156 


34 


16 


106 


sept.2004 


16 


E16 


55Blob 


1.716 


2.881 


151 


2 


14 


135 


sept.2004 


17 


Dl 


17Disk 


0.518 


0.173 


175 


12 


112 


51 


sept.2004 


18 


D2 


ISDisk 


0.693 


0.739 


203 


22 


12 


169 


sept.2004 


19 


D3 


52Disk 


0.715 


0.896 


207 


19 


17 


171 


sept.2004 


20 


D4 


51Disk 


0.879 


1.085 


190 


19 


22 


149 


sept.2004 


21 


D5 


04Disk 


0.862 


1.272 


204 


22 


19 


163 


sept.2003 


22 


SI 


sOl 


1.993 


-3.965 


74 


6 


32 


36 


sept.2002 


23 


S2 


s02 


1.740 


-3.524 


74 


3 


44 


27 


sept.2002 


24 


S6 


s06 


0.712 


-1.755 


89 


27 


44 


18 


sept.2002 


25 


S8 


s08 


0.189 


-0.945 


187 


69 


83 


35 


sept.2002 


26 


S24 


s24 


1.370 


-4.954 


131 


5 


35 


91 


sept.2003 


27 


S26 


s26 


0.593 


-2.809 


147 


5 


64 


78 


sept.2003 


28 


S27 


s27 


1.133 


-2.566 


145 


4 


77 


64 


sept.2003 


29 


D6 


lOlDiH 


-0.903 


-1.353 


279 


12 


257 


10 


sept.2004 


30 


D7 


102DiH 


-0.903 


-1.243 


270 


13 


233 


24 


sept.2004 


31 


D8 


104DiH 


-0.765 


-0.923 


273 


14 


218 


41 


sept.2004 


32 


D9 


105DiH 


-0.716 


-0.786 


273 


13 


244 


16 


sept.2004 


33 


DIO 


106DiH 


-0.693 


-0.605 


269 


15 


211 


43 


sept.2004 


34 


E22 


107ExH 


-1.397 


-1.643 


268 


5 


221 


42 


sept.2004 


35 


E23 


lOSExH 


-0.536 


-1.891 


259 


13 


185 


61 


sept.2004 


36 


E24 


109ExH 


-0.531 


-1.470 


251 


27 


174 


50 


sept.2004 


37 


E17 


llOHaS 


-2.177 


-2.402 


157 


27 


10 


120 


sept.2004 


38 


E18 


lllHaS 


-2.108 


-2.198 


123 


15 


10 


98 


sept.2004 


39 


E25 


123G1S 


3.080 


-3.161 


104 


3 


20 


81 


oct.2005 


40 


E26 


124G1S 


3.048 


-3.127 


109 


4 


21 


84 


sept.2005 


41 


D12 


131DiH 


0.653 


0.852 


289 


3 


39 


227 


oct.2005 


42 


D14 


134DiH 


0.867 


1.115 


279 


27 


25 


227 


oct.2005 


43 


D15 


135DiH 


1.189 


1.300 


213 


74 


15 


124 


oct.2005 


44 


E27 


137DiH 


1.338 


1.581 


262 


99 


21 


142 


oct.2005 


45 


E28 


138DiH 


1.356 


1.656 


265 


20 


26 


119 


oct.2005 


46 


E21 


148ExH 


-1.080 


-2.045 


213 


17 


116 


80 


sept.2005 


47 


E29 


150ExH 


1.206 


2.205 


229 


74 


38 


117 


sept.2005 


48 


E30 


151ExH 


1.695 


2.396 


195 


43 


16 


136 


sept.2005 


49 


E31 


152ExH 


1.679 


2.536 


209 


33 


18 


158 


sept.2005 


50 


S28 


153ExH 


0.1916 


-0.863 


170 


21 


84 


65 


sept.2005 


51 


D13 


165DiH 


0.747 


0.948 


265 


27 


16 


222 


oct.2005 


52 


Dll 


166DiH 


-0.6643 


-0.553 


240 


20 


132 


88 


oct.2005 


53 


E20 


167HaH 


-1.597 


-1.838 


209 


16 


128 


65 


oct.2005 


54 


E19 


168HaH 


-1.605 


-2.180 


170 


19 


35 


116 


oct.2005 


total 










9861 


1207 


4032 


4622 





abed 



"^Total number of stars with cross-correlation peaks greater than 0.05. 

''halo stars: all stars lagging the peak of the disk velocity by vi^g > 160km/s, or when no disk is 
detected in outer halo fields, then from the average disk lag, 25km/s. From this sample of 1207 RVs, 
the giant stream is windowed out, and only stars with velocity errors < 20km/s are used, giving a 
final Halo sample of 827 stars. 

"^Disk stars (encompassing both thin and thick components, as described in the text): stars rotating 
with our disk model with a velocity lag vi^g < 160km/s. 

''Stars with radial velocities falling within the Galactic window, defined as vjiei > — 160km/s. Note 
that this column is only for consistent book-keeping of the Halo sample across fields, but many stars 
from this window are used for disk analyses (see e.g., 105, and fig.|3. The actual Galactic fraction of 
the total sample is estimated at ~25%. 



